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The fluctuating envelope of the pulse-to-pulse radar echoes from a range 
cell consisting of a stationary target along with many independent, randomly 
moving scatterers is assumed to behave like a stationary Rayleigh process. 
In radar terminology this fluctuating or fading envelope of the pulse-to- 
pidse radar echoes is called signal plus clutter. The envelope of the pulse-to- 
pulse radar echoes may fade below some critical threshold level for a dura- 
tion such that the performance of the radar becomes unsatisfactory. Theoret- 
ical approximations for the probability densities of both the duration of fades 
and the interval between fades of the underlying Rayleigh process are pre- 
sented in graphs for various threshold levels and various signal-to-clutter 
power ratios. The corresponding exact residts are at present unknown. The 
results of this paper apply to all other fields of science and technology for 
which a stationary Rayleigh process characterizes a fading phenomenon. 

I. INTRODUCTION 

Consider a pulsed radar system "viewing" a range cell consisting of a 
stationary target along with many independent, randomly moving 
scatterers as shown in Fig. 1. Each received echo consists of the vector 
sum of all the elementary echoes originating from within the range cell. 
The contributions from the randomly moving scatterers arrive at the 
radar receiver with random phases. As a result, each received echo will 
consist of a steady signal, from the stationary target, plus a Gaussian 
perturbation. Accordingly, samples of the envelope of the pulse-to-pulse 
radar echoes can be considered as samples of an underlying Rayleigh 
process. Thus, the effect of the randomly moving scatterers is to cause 
the envelope of the pulse-to-pulse radar echoes to fluctuate or fade in an 
irregular manner. The envelope of the pulse-to-pulse radar echoes may 
fade below some critical threshold level for a duration such that the 
performance of the radar becomes unsatisfactory. 
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Fig. 1 — A model for studying the duration of fades associated with signal 
plus clutter. R{t,a) represents the envelope of the pulse-to-pulse radar echoes 
from the range cell. At the level R, <f> and <f> + represent the duration of a fade 
and the interval between fades, respectively. 

In radar terminology the fluctuating or fading envelope of these pulse- 
to-pulse radar echoes is called signal plus clutter. Classical discussions 
of signal plus clutter were given by H. Goldstein and A. J. F. Siegert 
and can be found in Refs. 1 and 2. A well-known example of signal plus 
clutter is the envelope of the pulse-to-pulse radar echoes from a target 
surrounded by a great deal of "chaff". Some other examples may be the 
envelope of the pulse-to-pulse radar echoes from the aurora, the iono- 
sphere, the sea, the ground, meterological precipitation, and a hyper- 
sonic object during reentry of the earths atmosphere. 

A natural assumption for studying the duration of fades associated 
with signal plus clutter is that the random process underlying the fading 
is a Rayleigh process. However, only a few theoretical results are avail- 
able concerning the duration of fades associated with Rayleigh processes. 
Thus, one is often unable to determine how well a Rayleigh process 
actually characterizes the duration of fades observed experimentally. 

The purpose of this paper is to present some additional theoretical 
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results which characterize approximately the duration of fades one would 
expect when the random process underlying the fading is indeed a Ray- 
leigh process. We shall assume that the envelope of the pulse-to-pulse 
radar echoes behaves like the Rayleigh process R(t,a) sketched in Fig. 1. 
R{t,a) represents the envelope of a stationary random process consist- 
ing of a sinusoidal signal of amplitude \/2a and frequency f plus a 
Gaussian process of unit variance having a narrowband power spectral 
density Wb(f — f ) which is symmetrical about f . We assume that the 
radar pulse repetition frequency is several times greater than the band- 
width of the Rayleigh process R(t,a) in order that an adequate number 
of radar echoes are used to form R(t,a). Also, we shall assume that the 
variance of the receiver noise is negligible in comparison with the vari- 
ance of the clutter. 

Using notation consistent with Refs. 3, 4, and 5 we shall present 
theoretical approximations for the following probability functions for 
arbitrary signal-to-clutter power ratio "a": 

(i) P ~(T,R,a), the probability density of the duration of a fade of 

the Rayleigh process below the level R. 
(ii) Pi(r,R,a), the probability density of the interval between fades 

of the Rayleigh process below the level R. 
(in) F ~(T,R,a), the probability that the duration of a fade of the 

Rayleigh process below the level R lasts longer than r. 
(iv) Fi(r,R,a), the probability that the interval between fades of the 

Rayleigh process below the level R lasts longer than t. 

The model considered in this paper also has application in the study 
of the duration of fades in radio transmission. In fact Rice 6 - 7 led the way 
by analyzing the duration of fades in radio transmission assuming that 
the underlying random process was R(t,0). 

II. INTEGRAL EQUATIONS AND EXPECTATIONS 

Let us define the following auxiliary probability functions for arbi- 
trary level R and arbitrary signal-to-clutter power ratio "a": 

CO Q~(T,R,a) dr, the conditional probability that an upward level- 
crossing occurs between t + r and t + t + dr given a downward 
level-crossing at t. 
(ii) [U(r,R,a) — Q(r,R,a)]dT, the conditional probability that an 
upward level-crossing occurs between t + t and t-\- r + dr 
given an upward level-crossing at t. 
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Explicit expressions for these auxiliary probability functions were 
presented in Ref. 5. 

Approximate theoretical results for P ~(t,R,o,) and Pi(r,R,a) are 
given by the following integral equations: 

P ~(T,R,a) = Q~(T,R,a) - P ~{t,R,<l) * [U(r,R,a) - Q(r,R,a)] (1) 

Pi(r,fl,a) = [U(T,R,a) - Q(r,R,a)] 

(2) 
- Pi(r,R,a) * [U(r,R,a) - Q(r,R,a)] 

where * denotes the convolution operator, that is, 

f*g = f f(t)g(r - t)dt. 

•'—00 

Equations (1) and (2) were derived in Ref. 3 by applying McFadden's 8 
"quasi-independence" idea to the Rayleigh process R(t,a). Also, by 
definition we have that 



F, 



' '(T,R,a) = 1 - / Pa (T,R,a)dr (3) 



and 



Fi(r^,a) = 1 - f Pi(r,/2,a)rfr. (4) 



The exact expectations E ~(t,R,ci) and Ei(t,R,cl) associated with the 
respective densities P ~(T,R,a) and Pi(r,R,a) can be computed from 
the following equations: 

, , _ PrjR(t,a) < R] 
K {T,R,a) = ^ 

00 (5) 

£ (R 2 "/2"nl) 1 F 1 (n+ ±; 2n + 1; -2R\/2a) 



V(/3/2t) RiFiih 1; -2RV2a) 

El (r,R,a) = 1 = exp [(R - V2aYj2\ (fi) 

* VW2T) RiFt (fc 1; -2R^2a) 
where 

= (2tt) 2 f W b (f - /.)(/ - fof df 

Wi(f — f ) = narrowband power spectral density of the Gaussian 
process involved in the definition of R(t,a) 
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average signal power 
average clutter power 

iFi(a;0;x) = the confluent hypergeometric function 
_ -, , a , a(a + 1) X 

~ ^ ^ 0(0 + D2! 1 " '"• 

Pr { } = probability of the event inside the brace 

N R = average number of upward (or downward) crossings of 
the level R per second. 

Equations (5) and (6) were developed in Ref. 4, and they follow directly 
from some well-known results reported by Rice and Bennett. Each 
iFi function appearing in (5) and (6) can be expressed in terms of a 
Bessel function of imaginary argument. 

Thus, with the aid of a digital computer one can compute theoretical 
approximations for the probability functions of interest in this paper 
along with the exact theoretical expectations given by (5) and (6). 

III. RESULTS FOR A GAUSSIAN AUTOCORRELATION FUNCTION 

In order to define the Rayleigh process R(t,a) underlying the fading 
phenomenon we need to specify both Wb(f — /©) and the signal-to- 
clutter ratio "a". The normalized autocorrelation function ra(r) asso- 
ciated with W b (f — }o) is given by 

m ( T ) = f W b (f - fo) cos 2a-(/ - /„)«*/. (7) 

Thus, ra(T), rather than W b (f — f ), can be used to define the Rayleigh 
process R(t,a) underlying the fading phenomenon. Notice that /3 appear- 
ing in (5) and (6) is merely — m" (0). The primes denote differentiations 
with respect to r. 

Ref. I points out that it is convenient to measure the normalized 
autocorrelation function of the fluctuating low frequency power P(t) = 
R 2 (t,a) and denotes this normalized autocorrelation function by p(P,r). 
In explicit terms p(P,t) is denned as 



p(P,r) = 



where 



E{[P(t + r) - EP(t)][P(t) - EP(t)}\ 

Var P(t) 

(8) 
_ EP(t + r)P(t) - E 2 P(t) 

Var P(t) 
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E = Expectation 
Var = Variance. 
Refs. 1 and 2 relate m (t) and p (P,t) as follows: 

m(r) = Va 2 + (1 + 2a) P (P,r) - a. 



(9) 



Ref. 1 also points out that the appropriate value of "a" can be estimated 
by measuring the probability density of P(t) and comparing the result 
with the theoretical probability density of P(t). Thus, (9) indicates 
that the Rayleigh process R (t,a) underlying the fading phenomenon can 
also be defined from measurements of the normalized autocorrelation 
function p(P,r) and the value of "a." 

For purposes of computation we shall take Wb(J — fo) and m(r) as 




4 5 6 

Ub=2-7rCTbT 



Fig. 2 — P<r{v.b , ko , a) is the probability density of the duration of a fade of 
the Rayleigh process below the normalized level ko . The autocorrelation function 
of the Gaussian process involved in the definition of the Rayleigh process is m(r) 
and the signal-to-clutter power ratio equals "a." 
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Fig. 3 — Pi(ub , ko , a) is the probability density of the interval between fades 
of the Rayleigh process below the normalized level ko ■ The autocorrelation func- 
tion of the Gaussian process involved in the definition of the Rayleigh process is 
m(r) and the signal-to-clutter power ratio equals "a." 



follows: 



and 



WW - /.) - 



Ob 



V'2w 



exp 



-(/-/. 



2<T6 2 



q 



»i(t) = exp 



xp 



— (27T(T6t)' 



(10) 



(11) 



This particular choice tends to characterize the radar clutter fluctua- 
tions observed experimentally. 1 ' 2 From (11) we see that it is convenient 
to define normalized time as Ub = 2wobT. 

For the experimenter it is convenient to normalize the threshold level 
with respect to the average value, ER (t,a), of the Rayleigh process. We 
shall consider three such normalized levels k 



R 



-*. = !, 
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T ' V2T ' 



ER(t,a) 

The expectation ER (t,a) was derived by Rice 9 and is given by 



(12) 
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ER(t,a) = j/ftfi^-ijl;-*). (13) 

When a = we have that R = V (71-/2), 1, $, These latter two values 
of R were also considered by Rice for the case a = 0. 

Figs. 2 through 10 present the computed results for a = 0, 1, 4, and 
k = 1, -\/2/t, l/\/2ir. The numerical evaluation of Q~(r,R,a) and 
U(r,R,a) — Q(r,R,a) was carried out by using Simpson's rule. Integral 
equations (1) and (2) were solved numerically by using the trapezoidal 
rule. All results are plotted with respect to normalized time Ub . The cor- 
responding experimental results for k = 1 and a = 0, 1, 4 were presented 
in Ref . 4, and they agree well with the approximate theoretical results 
presented in Figs. 2, 3, and 4. 
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Fig. 4 — F ~(ub , k , a) is the probability that the duration of a fade of the 
Rayleigh process below the normalized level & lasts longer than Ub . Fi(ub , k , a) 
is the probability that the interval between fades of the Rayleigh process below 
the normalized level &o lasts longer than Ub . 
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Fig. 5 — Pd~(ub , k , a) is the probability density of the duration of a fade of 
the Rayleigh process below the normalized level &o . The autocorrelation function 
of the Gaussian process involved in the definition of the Rayleigh process is m(r) 
and the signal-to-clutter power ratio equals "a". 
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Ub=27TCTbT 

Fig. 6 — Pi (ub , ko , a) is the probability density of the interval between fades 
of the Rayleigh process below the normalized level ko . The autocorrelation func- 
tion of the Gaussian process involved in the definition of the Rayleigh process is 
to(t) and the signal-to-clutter power ratio equals "a". 
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Fig. 7 — Fo~(ut ,k ,a) is the probability that the duration of a fade of the 
Rayleigh process below the normalized level k lasts longer than Ub ■ Fi(ub , k ,a) 
is the probability that the interval between fades of the Rayleigh process below 
the normalized level ko lasts longer than Ub • 



For deep fades and large signal-to-clutter power ratio "a", one would 
expect P ~(r,R,a) to approach a Rayleigh probability density. For as 
"a" gets large the Rayleigh process R(t,a) tends to behave much like a 
Gaussian process, see (3.6) of Rice, 9 and the durations of deep fades of 
Gaussian processes are known to be characterized by a Rayleigh proba- 
bility density. 7,10 Figs. 8 and 10 show that this is, approximately, the 
case when k = l/y/2r and a = 4. Thus, for A; ^ 1/-\/2tt and a ^ 4 we 
have the following approximate results: 



Pr(r ^)= ' (^1, 



7T 



m 



(14) 



and 
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Fa 



!r^> -«p[-|(£)"]. (15) 



The value of E appearing in (14) and (15) is given by (5) with 
R = k„ER{t,a). 

Equations (14) and (15) are useful approximations when k is small 
and "a" is large for an arbitrary normalized autocorrelation function 
w(t) such that m" (0 + ) = 0, although we have been treating the re- 
strictive Gaussian autocorrelation function defined by (11). The condi- 
tion m (0 + ) = leads to Q~ (0 + ,R,a) = 0, and thus the approximation 
given by (14) is exact at t = + . As a partial check on this gener- 
alization we also verified that (14) and (15) begin to be useful approxi- 
mations when k = l/x^rr, a = 4 for the normalized autocorrelation 
functions 
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Fig. 8 — Po~(ub , ko , a) is the probability density of the duration of a fade of 
the Rayleigh process below the normalized level ko . The autocorrelation func- 
tion of the Gaussian process involved in the definition of the Rayleigh process is 
?n(r) and the signal-to-clutter power ratio equals "a." 
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and 

m ( T ) = ^l + cd2 | T | + i^]exp(-co 2 |r|). (17) 

Equation (16) corresponds to an ideal bandpass power spectral density 
Wtif - fo) given by 

((2/c)- 1 for fo-fe£f£fo + fe 

W b (f-fo) = \ < 18 ) 

\0 otherwise. 

Equation (17) corresponds to a power spectral density W b (/ — /<>) 
given by 

8/(3tt/ 2 ) 



where 



W b (f - f ) = 



[>*mr 



(19) 



cu2 = 2irj2 . 
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Fig. 9 — Pi (ub ,k ,a) is the probability density of the interval between fades 
of the Rayleigh process below the normalized level h . The autocorrelation func- 
tion of the Gaussian process involved in the definition of the Rayleigh process is 
7ji(t) and the signal-to-clutter power ratio equals "a". 
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Fig. 10 — F ~(«6 , ko , a) is the probability that the duration of a fade of the 
Rayleigh process below the normalized level kt> lasts longer than Ub . Fi(ub , ko , a) 
is the probability that the interval between fades of the Rayleigh process below 
the normalized level k lasts longer than Ub ■ 

For a given m(r) with ra'"(0 + ) = along with a ^ 4, k = 1, the 
duration of fades and the interval between fades of the Rayleigh process 
R(t,a) behave as if they were generated at the mean value level of a 
Gaussian process having a normalized autocorrelation function of mir). 
For example, compare P a ~(ub , k„ , 4) of Fig. 2 with the experimental 
points plotted in Fig. 2 of Ref. 3. Also compare Pi(w& , k , 4) of Fig. 3 
with the experimental points plotted in Fig. 3 of Ref. 3. 



IV. CONCLUSIONS 

Assuming that the random process underlying a fading phenomenon 
is a stationary Rayleigh process, one can compute useful theoretical 
approximations for the probability functions which characterize the 
duration of fades and the interval between fades. The corresponding 
exact results are at present unknown. 



1298 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1966 

For deep fades and large signal-to-clutter power ratio the duration of 
fades is characterized, approximately, by a Rayleigh probability density. 

For large signal-to-clutter power ratio the duration of fades and the 
interval between fades of the Rayleigh process below the mean value 
level behave as if they were generated at the mean value level of a 
certain Gaussian process. 

The results of this paper apply to all fields of science and technology 
for which a stationary Rayleigh process characterizes a fading phe- 
nomenon. 
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